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ABSTRACT 



Three previously unpublished methods for ertpirically detentdning 
deep-water sv;ell direction were examined in this vrork: (1) triangulation 

using two widely separated wave sensors, (2) intersection using weather 
maps and a single wave sensor, and (3) swell point-source estimation frcm 
weather maps. The primary objective of each method Wcis to identify a 
single point source of the swell train produced in an approaching cyclonic 
storm. Method (1) was not adequately tested, but results as applied to 
swell fron one storm were favorable. Methods (2) and (3) , applied to 
five selected North Pacific storms, gave close agreement on the swell 
origin time (15 itiinutes to 6 hours) and source point location (27 to 362 
nautical miles) , and on swell arrival direction at Monterey, California 
(0.1 to 5.3 degrees) . Methods (1) and (2) give deep-water directions for 
swell trains that have already arrived, while Method (3) is a prediction 
method. 
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I. INTPCBUCTiaJ 



Deep-water swell directicm information is needed for most operational 
coastal engineering work. Swell direction data, used by coastal engineers 
for such purposes as the design of structures and the coi^tation of 
littoral dirft, can be found in the form of statistical corpilations 
derived from wave-hindcasting methods. However, the deep-water direction 
increments used in such corpilations are 45 degrees and 22.5 degrees. 

These direction intervals are crude vdien the sensitive dependence of 
refraction on deep-water wave direction is taken into account in conputing 
shoal-water wave data for a selected coastal site. Swell direction, 
generally in the form of a mean direction, can be obtained from wave fore- 
casting or hindcasting from weather maps , from multiple sensor arrays 
placed in shoal water, and fron more exotic procedures such as sea-glitter 
measurements and laser profiling from aircraft. The latter methods are 
not used to obtained operational wave-direction information, and have the 
disadvantage of being very costly. 

In this study three previously undescribed methods for obtaining deep- 
water swell direction were investigated using recorded wave data fron 
coastal sensors located at Monterey and San Clemente Island, California 
and synoptic weather maps. The three nethods were tested on five selected 
storms occurring from October through Decanber, 1971 \drich were tracked 
across the North Pacific. The swell direction from each storm was 
determined for deep water directly off Monterey. 

The three methods are all based on the concept that the doninant swell 
observed at a distant location is generated at a point source in space and 
time. They all involve identification of the point source. Once this 
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swell source is determined, the swell are then follcwed at group velocity 
to the coastal station of interest along a great-circle path, and the 
azimuth of the path at the coastal site is taken to represent the deep- 
water arrival dii'ection for the swell train. Ihe single direction that 
is obtained is considered here to be the direction of the dominant swell 
arriving at the station. It is recognized that real swell display both 
a directional and frequency spectrum; however, it is also evident that 
the size of ocean storms is normally such that for most storms at distances 
on the order of 1000 nautical miles or greater the directional bandwidth 
is narrow and sharply peaked. Accordingly, the determined azimuth is 
considered to represent the dominant band of the arriving directional 
spectrum with good precision. 

In the present study, a point source is thought of as that region 
of the sea surface (small ccnpared to the travel distance) and that time 
interval (small corrpared to the travel time) in which a prcminant energy 
itaximum is generated by the storm and from which the energy is then 
propagated toward the distant station of interest. It is assumed in 
this work that the swell observed at a distant station then represents 
the dispersive arrival of a burst of energy irtparted to the sea at the 
source location and time. lb identify this point source the three 
approaches formulated and tested were: 

1. Triangulation using swell source distances derived frcm analysis 
of the wave records frcm twD coastal sensors at widely separated 
locations . 
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2. Use of a single coastal wave sensor and synoptic weather maps. 

Ihe sensor gives the swell origin distance and time, and the weather maps 
give a probable path of the point of maximum wave generation toward the 
coastal station fron map to irap during the life of the storm. The inter- 
section of the distance arc and the path of maximum wave generation 
represents the single source point for the swell. 

3 . Determination of the swell origin point fron weather maps , by 
examining the surface pressure fields, in order to find that set of wind 
conditions which produces the maximum wave energy directed toward the 
coastal station. 

The North Pacific storms and swell trains, occurring in late 1971, 
examined in this paper, are indicated as follcxvs: 



STOKVSWELL TRAIN 
NUMBER 


ST0H4 

OCCURENCE 


SWELL ARRIVAL 
AT MONTEREY 


1 


25-26 October 


30 Oct - 1 Nbv 


2 


31 Oct - 5 Nov 


7-9 Noverber 


3 


19-24 November 


25-26 November 


4 


22-26 November 


28-29 November 


4A 


26-29 November 


29 Nov - 1 Dec 


5 


11-14 Deconber 


16-18 December 



Storms 1, 2, 3, and 5 were unccstplicated storms that developed and 
traveled across the North Pacific much in the manner expected of winter 
storms. Storm 4 presented an analysis problem. After several days of 
development it weakened considerably and appeaured ready to die out. 
However, the center of the low shifted its position about 420 nautical 
miles to the southwest over a 12-hour map interval, and the storm 
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re-intensified. After about 24 hours the storm then moved to the east 
as normally expected. For this reason, and because spectral analysis of 
the waves at Monterey indicated two arriving swell trains, it was decided 
to analyze this particular event as two separate storms, each with a 
unique source point. The earlier portion of the storm was accordingly 
labelled Storm 4, and the' latter Storm 4A. 
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II. SWELL GENERATION AND PROPAGATION CONSIDERATIONS 



The five cases studied all deal vd.th the moving fetch problem of 
swell produced in a moving cyclonic storm where the swell is propagated 
toward a coastal station located approximately in the direction of the 
sotrm's advance. In this situation, it is usual for the wave-generating 
area to travel for a number of hours in the same direction and at the 
same speed as the group velocily^ of some of the frequency ccnponents pro- 
duced in the sea. Accordingly, the effective dxoration of the wind over 
the waves is long, the maximum amount of energy the prevailing wind can 
transfer to the sea is transferred, and the sea approaches or reaches a 
fully arisen condition. 

The maximum energy put into the sea during the passage of a storm 
across the North Pacific, in the case of waves traveling in the direction 
of the storm, is determined by two factors: (1) the wind speed, and 

(2) the relative speed of the wave generating area with respect to the 
speed of the wave energy propagation. It can be expected that in storms 
where the difference between the speed of the storm and waves is always 
large, the maximum energy will be put into the waves about the time of 
maximum winds. On the other hand, the tirre of rraximum energy input into 
the sea in a moving storm can be expected to occur viien the storm speed 
equals the group velocity of the waves of the peak frequency that can be 
generated by the wind then blowing, even through the winds may have been 
stronger at a time viien the relative velocity of storm and waves was very 
different. 
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During its passage across the North Pacific, the speed of the wave- 
generating area is usually roughly the same as the group velocity of the 
longer waves that can be generated by the winds present. It thios appears 
that the sea generated in the direction of storm travel must reach a 
fully arisen condition, either when the wind in the storm reaches its 
peak mean velocity or when the relative speed of the storm and waves is 
coincident. That time is then the time when the energy in the sea is at 
a maximum. 

One of the fundamental relationships used in this study, both in the 
analysis of wave records to determine source distance and origin time of 
arriving swell trains and to predict fron weather maps the arrival tiroes 
and frequencies in a swell train generated in a cyclonic storm, was the 
frequency- time (f-t) graph introduced by Monk, et al. (1963) . 

Using spectral analysis of successive wave records, these investi- 
gators determined the energy density at selected frequencies and constructed 
frctn this information a frequency-tirte graph showing the spectral energy 
density distribution by means of contours for the period of time covered 
1:^ the wave records. These graphs show that a swell train arriving from a 
storm is revealed as a prcrninant sloping ridge line in the spectral density 
topography. For swell frcm storms at moderate and great distances these 
ridge lines are linear. 

In the context of linear wave theory, a straight f-t line inplies a 
point source for the swell train in both space and time, and provides a 
simple means for determining the apparent generating point location and 
origin time. In practical application a point source implies that the 
linear dimensions of the wave generating area are small compared to swell 
travel distance and that the source duration is small ccmpared to travel 
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time. In the present study, swell fran soiorce distances as close to the 
wave recording station as 1970 nautical miles were treated. The swell 
signatures on the f-t graphs constructed using recorded wave data were 
linear for these distances, and the source infomation derived for the 
swell train was consistant with the weather maps. It may be noted here 
that Snodgrass, et al. (1966) dealt with swell sources at distances as 
small as 1320 nautical miles. 

The distance and origin time of the swell source with respect to the 
wave sensor can be obtained from the f-t curve fitted to the swell data 
in the following way. If the energy associated with each frequency com- 
ponent, f , generated in the sea is propagated at its own group velocity, 

C , over a distance, D, from the generating point to the observation 
point in travel time, t-t^, then 

Cg - D/(t-t^) 

but 

Cg = g/4irf 

from linear wave theory, so that 

g/4Trf = D/ (t-t^) . 

The slope of the f-t line is then 

df/dt = g/4TrD. 

Therefore, the apparent distance fixm the estimated source point is given by 

D = 1.515/ (df/dt) nautical miles (1) 

where f is in Hertz and t is in hours. The estimated origin time, t^, 

is given by the intercept of the f-t curve with the zero frequency line. 
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Equation (1) assumes that the distance, D, between the source and 
observation point is constant. If this were not the case, there would 
be an additional term in Equation (1) to account for the change in slope 
of the f-t curve as the storm moves across the ocean. The assumption of 
Monk, et al. (1963) that the f-t curve is linear is interpreted in this 
thesis to mean that the energy is iit^ulsively transmitted to the waves 
at the space-time point-source; therefore, the distance of swell propa- 
gation to a distant station may be considered constant for this burst of 
energy. 

The specific applications of the f-t graph in this study will be 
shewn later. 
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III. DETERMINAriaN OF EEEP-WMER SWEUD DIPECTIQN 



A. METHOD ONE: TRIANGULATICN USING IWD SENSORS 

Ihis method involved triangulation to determine the swell source point 
losing data obtained from spectrally analyzed wave records of a given swell 
train recoorded at two widely separated coastal wave sensors. The source 
point was located by striking arcs of great-circle distance from each 
sensor site, the radii of vdiich equal the swell soiorce distances as 
derived fixm the f-t curve for each station. 

The original plan was to use the wave gages at San Clemente Island, 
California and Newport, Oregon, vAiich are about 700 nautical miles apart; 
however, failure of the Newport sensor cable precluded its use. The 
Monterey gage was chosen as an alternate sensor. This choice reduced the 
wave-gage separation to about 300 nautical miles. It also placed the 
Monterey and San Clemente Island sensors roughly on a line with the swell 
source point to be located, thereby reducing the intersection angle 
beto\7een the two radii to a small value. Because of this unsatisfactory 
triangulation geonetry, this method was applied to the swell from one 
storm only. Storm 3 was chosen because of its high winds and unccarplicated 
character. 

1. Wave Recording and Spectral Analysis 

The arriving waves at Monterey were recorded by a Marine Advisers 
Model A-2 pressure sensor located in 28 feet (MWL) of water off Del Monte 
Beach. Waves arriving at San Clemente Island were recorded on a pressure 
sensor located near the northwestern tip of the island in 40 feet (MWL) 
of water. This sensor is naintained and operated by the Fleet Weather 
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Facility, San Diego, California. The waves at both gages were recorded 
in analog form on a paper strip chart. The records consisted of a 
twenty-minute fast trace every six hours recorded at a chart speed of 
one inch in 30 seconds, separated by a slow trace. 

The consecutive twenty-minute fast traces for the duration of the 
swell train at both stations were first converted from analog form to 
digital data using a Calma Ccmpany digitizer made available by the Fleet 
Numerical Weather Central (FNWC) , Monterey, California. This taped 
information was transferred to punched cards using the program CONVERT 
(Lynch, 1970) designed for use in the FNW3 CDC 6500 ccmputer. 

A frequency spectrum analysis was then performed on the digitized 
wave record using the Blackrtan-Tukey (1958) method. As evidenced by 
the spectra produced, the highest wave frequency measured was about 0.20 
Hz; accordingly, the number of degrees of freedom was 14.5 and the 90 
percent confidence limits were 0.60 and 2.10 of the energy density 
estimates for the width of spectral bands chosen. The output of the 
analysis was wave frequency (Hertz) and the corresponding energy density 
(feet^-second) in both tabular and graphical forms. The energy values 
were ccmputed for frequencies from 0.0 to approximately 0.24 Hz, in 
0.00416 Hz increments. 

2. Determination of Swell Source Distance and Origin Time 

The spectral energy density at six-hourly intervals during the 
period of swell arrival was plotted on a frequency-time graph for each 
station. An f-t line was then fitted to the spectral energy peaks. Fran 
the slope of the f-t curve the distance to the point of origin was 
determined, and fron its intercept with the zero frequency axis the time 
of origin was determined, as explained above. The frequency-time graph of 
spectral energy density for Swell Train 4 is shown in Figure 1. 
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It is possible for the frequency of maximum energy density in 
swell to shift across the spectrum due to shoaling effects as the waves 
travel frctn deep to shallow water, Thanpson (1970) investigated this 
phoiamenon and determined that the magnitude of a frequency shift of 
this sort is negligible for swell passing over a narrow coastal shelf. 

3. Great-Circle Distance and Trajectory Overlays 

In order to quickly and easily measure great-circle distances on 
the weather maps, to determine the trajectory frcm the swell source to 
the station, and to measure the arrival direction of swell at the station, 
a caiputer program was written using a Naval [Research Laboratory sub- 
routine titled TO-NE?L-GCDIST (Chang, 1969) . The program produced a 
tabulation of great-circle distances and the azimuths of these great- 
circles frcm Monterey and San Clemente Island to integer latitude- 
longitude intersections (one-degree increments) in the North Pacific 
Ocean. Using this table, two very convenient transparent overlays were 
drawn for each station for use on the PTSIWC strip chart of the North 
Pacific Ocean (1:30 million-scale polar stereographic projection) . The 
first of these overlays contains a family of arcs of great-circle dis- 
tance in 100 nautical-mile increments, vhile the second shows great-circle 
paths drawn for five-degree increments of arrival direction at the 
station. Distance overlays were drawn for Monterey and San Clemente 
Island. These overlays and the associated table were used extensively 
throughout this investigation. 
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4. /^plication of Methcxi One 



As mentioned above, this method was applied only to Swell Train 
3, viiich was observed to arrive at the two coastal stations on 25-26 
NOvaiiber , 1971 . 

The source distances and origin times obtained from the f-t 
curves drawn from the swell data at the two stations are as follovs: 



Station 


Source Distance 


Origin Time 


hm 


2330 naut mi 


032/21 Nov 71 


SCI 


2640 naut mi 


002/21 Nov 71 



Triangulation was performed by constructing arcs frcm both gages 
at their respective source distances on a base map using the distance 
overlays. These arcs are shown in Figure 2. 

Triangulation cannot be performed directly using these arc dis- 
tances because the origin times are not coincident. In order to resolve 
this problem one of the sets of data was accepted as a reference, namely 
that for the Monterey gage. The data frcm the San Clemente Island gage 
were then adjusted as follows. Ihree hours were added to the origin time 
for the swell arriving at San Claiiente Island to make it coincident with 
that at Monterey. This time adjustment represents a closer swell source 
for the San Clemente Island site than was obtained frcm the f-t graph in 
an amount given by the distance waves would travel frcm the source in 
three hours. Because many frequencies are generated by the storm, each 
of which travels at a different group velocity, the frequency of maxirautn 
energy density in the swell arriving at San Clemente Island = 0.07 Hz) 
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was used to give a representative travel distance. Thus, the travel 
distance for the three-hour incranent. At, is: 

D = C At = = 66 nautical miles. 

g f 

The adjusted swell source distance from San Clanente Island is then 
2640 - 66 = 2574 nautical miles. This adjusted distance is shown in 
Figure 2 by the dashed line. 

Because the gages at the two stations are nearly in line with the 
swell source, as is evident frcam the figure, a satisfactory intersection 
angle between the two arcs could not be obtained. However, the generally 
good source distcince agreement between the dashed arc drawn fron the 
San Clo:iiente Island gage site and the solid arc drawn fron the Monterey 
site is apparent, and was found to agree with a storm location on the 
weather maps, as discussed later. 

B. METHOD TWD; TTUTySIGUIATION USING WAVE RECOBDS AND WEATHER MAPS 
1. Procedure 

The second method for corrputing the deep-water swell direction 
involves detemdnation of the point of swell origin \'d-thin the storm 
using the intersection of a swell source distance arc as determined by 
analysis of the wave records from a coastal wave sensor with the path of 
the peak winds (defined below) directed toward the sensor site as obtained 
fron analysis of a series of synoptic weather maps. The swell arrival 
direction is then detennined fron the great-circle trajectory drawn from 
the intersection point to a coastal site at or in the vicinity of the 
wave gage. 
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Ihe distance from the wave gage to the origin point of the swell 
is determined from the f-t curve for the swell train by the procedure 
described under Method One. This distance is then struck as an arc from 
the sensor site on a blank weather map used for plotting of the data. 

The maximum wave energy generated by the storm in the direction 
of the wave-gage site is considered to be that generated along the peak 
wind path. Ihe path of maximum sxirface wind speed directed tcv/ard the 
wave gage site was obtained as follows; 

For a. given storm, the path of peak winds directed toward the 
wave-gage site was obtained by plotting the peak wind point derived from 
each surface pressure analysis onto the blank data itap, and then 
connecting these points of maximum surface winds directed toward the 
station as is shown in Figure 3 . The point of effective maximum winds 
on an individual weather map was obtained by overlaying on the itap the 
great-circle path terrplate for the Monterey gage, then by using a second 
transparent teirplate indicating a fifteen degree cross-isobar surface 
wind flow, locate the axis of the storm along which the surface winds 
blew 15 degrees across the isobars and directly toward Monterey. A 
constant fifteen-degree cross-isobar flow angle for the surface wind was 
chosen as an average value based on cormion usage. A foxor-millibar isobar 
spacing, chosen to give the maximum wind speed, was then measured along 
this axis and conveirted to the geostrophic wind velocity, V . The point 
of maximum vnnd velocity on the axis was considered to represent the 
source of maodmum wave energy being generated on that map in the direction 
of Monterey. The point of maximum wind speed on each weather map 
obtained in this way was then connected to produce the path of peak winds 
for the storm. 
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In cxJitplicated wind fields, wave energy that is directed toward 
the station fron locations off the peak wind axis may represent the 
maximvim energy contributed from a given weather map. This would lead 
to distances and origin times somewhat different from that obtained from 
this procedure, 

2. Example of Method IWo Application 

Storm 5 is presented as an example of the procedures used in 
Method 'IV;o, This particular storm was first noted on the weather maps 
of 11 December 1971, and the first significant swell arrival at Monterey 
from wave record analysis was can 16 December 1971 

Figure 4 shows the FNW2 surface pressure map for OOZ 13 Dec 71 
analyzed to shew the axis of the surface winds in the storm that are 
blowing 15° across the isobars and directly toward Monterey. The circle 
represents the point of maximum winds on this axis, and the cross the 
center of the low-pressure system. The dashed line represents that por- 
tion of the axis which lies over and to the north of the Aleutian Island 
chain, and along which waves generated will not reach Monterey. The wind 
profile along the axis shown in Figure 4, which is typical of the storms 
studied, is shown in Figure 5. 

Figure 3 shows the axes of winds blowing toward Monterey as 
analyzed from eight successive 12-hourly surface weather maps. The 
circled points on the axes represent the points of maximum winds on each 
map, and are connected to give the path of maximum winds. On this map, 
the great-circle distance arc (heavy dashed line) is drawn at a distance 
of 2177 nautical miles, as computed from the f-t diagram for Swell Train 
Five. The intersection of the distance arc and the path of maximum sur- 
face winds (large circle) marks the presumed swell source point. The 
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Figiare 3: ' ANALYSIS OF STOFM 5 USING METHOD TVTO 
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Figure 4; FNWC (XiMPUTER-DRAWN SURFACE PRESSURE ANALYSIS FOR OOZ 13 DEC 71 SHCWING AXIS OF PEAK WINDS. 
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lATITUDE (DEGREES NORTH) 

Figure 5: WIND PROFILE ALONG AXIS OF PEAK WINDS 

ON OOZ 13 DEC 71 SURFACE PRESSURE MAP. 
(Storm 5) 
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Table I: RESULTS OF METBDD TWO 



L 

O (D 
(D 5-1 
U ( 1 ) 
•H -P 

o c 



•H 



o 



ro 


m 


in 


O 


(N 


ro 


• 


• 


• 


• 


• 


• 




00 


VD 


1 — 1 


ro 


o 


o 


o 


00 


o 


00 


1 — 1 


ro 


fO 


CN 


m 


CN 


ro 



4J 

•S« 

S3 

(D C7> 

p C 

I' 

M-l 

o 






W 




12 


W 




:s 


o 


in 


o 


in 


o 


o 


• 


• 


• 


• 


• 


• 




o 


1 — 1 


<Tl 


00 


a\ 


VD 




r- 


VD 


VD 


VD 


1 — ! 


I — ! 


1 — I 


1 — 1 


1 — 1 


•H 



SI 




IS 


IS 


IS 


S 


IS 


•d.ti 


o 


o 


in 


in 


o 


o 


• 


# 


« 


# 


• 


• 


rd -P 


as 


rH 






VD 


CN 


0 d 




in 


ro 




ro 


in 



*n 




1 — 1 


1 — 1 


iH 


iH 


rH 


1 — 1 












r- 




r- 


q 


1 — 1 £ 


-P 


> 


> 


> 




0 


•H 


1 — I *H 


0 


o 


o 


JQ 


Q 






0 ) E-» 


<D 


s 


IS 


IS 




Q 


•H 
















M 


(55 q 


in 




i-H 




VD 


m 


0 


•H 


CN 




CN 


CN 


CN 


rH 






N 


N 


N 


[S3 


N 


t<3 


e i! 


O 


O 


O 


O 


in 


in 




A4 o 


o 


o 


m 


o 


1 — 1 


rH 






in 


VD 


m 


o 


o 


1 — 1 


(D 




! — ! 


O 


o 


o 


CN 


o 



I 

-p 

CO 

8 



0) ^ 

Is 

CO 



4J 



o 


o 


o 


CN 


o 


r- 


r-H 


CN 


(T) 


00 


o 




O 


CN 


cn 


o 


CN 


1 — 1 


ro 


CN 


CN 


cn 


CN 


CN 



<N 



ro 









lO 



28 



great-circle propagation path of the swell of peak energy from the source 
point to Monterey is shewn by the heavy solid line. The deep-water 
arrival direction at Monterey is from an aziinuth of 310 degrees (true) . 

3. Results of Ifethod TWo 

Ihe results of this method of locating the point source and 
determining the swell arrival direction for all five storms are presented 
in Table I. Mditionally, Figure 6 shews the f-t curves derived for each 
of the five storms studied. 

C. METHOD THREE: PREDICTION FRCM WEATHER MAPS 

Method Three is an arpirical approach to identification of those 
events in a storm history which will allcw a forecaster to predict 
accurately the desired swell arrival information. This method uses 
surface pressure analysis weather neps to estimate the time and point of 
swell origin, and hence the arrival direction of the peak swell energy 
at a selected coastal station. The peak svell is that swell which is 
considered to be produced at the space-time point source by the impulsive 
introduction of energy into the sea. 

Mditionally, a useful extension of the data <±)tained using this 
method allows one to predict the time of arrival of the peak swell at 
the station. 

1. Time of Origin and Source Point of the Peak Swell 

The principal task in the analysis is to identify the time of 
origin and associated source point of the swell of maximum energy produced 
by the storm that is directed toward the coastal station. The individual 
weather maps shewing the storm of interest are examined, and their peak 
wind points are plotted on a blank base map, as outlined in Method IWo 
(illustrated in Figure 3) . The connected path of points of peak wind 
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velcx:ity then represents the Icxus of possible peak-swell origin points 
for that storm. To determine the swell origin point on the path of peak 
winds, one must next estimate the ircment in the storm history when the 
maximum amount of energy was imparted to the sea. In this work it was 
assumed that, in the absence of obvious overriding conditions such as 
higher velocity surface winds at an earlier time with little or no storm 
movement, the time of origin of the peak swell (and of a fully arisen sea 
condition) is that time vhen the group velocity, C^, of that frequency 
corrponent of maximum energy that can be produced in a fully arisen sea 
by the winds present equals the velocity ccmponent, V^, of the peak 
wind toward Monterey. 

A graph showing the history of and with time was prepared 
for each storm, as shewn in Figure 7 for Storm 5.. That monent vhen 
Cg = is considered to be the source time. It may be seen in Figure 7 
that this condition is met on both 11 and 13 December. The 13 December 
time was considered to be the most probable origin time because of the 
slightly stronger surface winds and closer approach of the storm to 
Monterey at that time. Once the morent in the storm's history vhen 
and were equal was determined, the location of the point source was 
established on the path of possible points on the base map by interpolation. 

Once the swell source has been chosen, the deep-water arrival 
direction of the swell is found by use of either the transparent great- 
circle path overlay or the associated table of values. 

The procedures for obtaining and are outlined as follows: 

(1) The velocity corrponent, V^, of the source point toward the obser- 
vation station is obtained by dividing the difference in the great-circle 
distance from the peak wind points on successive weather maps to Monterey 
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by the time interval between the maps. (2) The group velocity, C , of 
the ccnponent of maximum energy is obtained fron the geostrophic wind 
velocity. The geostixjphic wind velocity, V , for the peak wind point is 
obtained fron the pressure gradient on each surface pressure map. This, 
in turn, is converted to the surface wind velocity, V^, using an empiri- 
cally derived ratio of V /V =0.66. This constant is for conversion of 

s g 

geostrojluc wind velocities obtained using a four-millibar isobar spacing 

(discussed under Method TVro) . Assuming that fully arisen sea conditions 

exist on all maps, the frequency on maximum energy, f , produced by 

the surface wind on each map was obtained from Figure 8 (the construction 

of this figure is based on the Pierson-Moskowitz (1964) spectrum) . The 

group velocity of the frequency coirponent of maximum energy is then found 

from the propagation equation, C = 1.515/f (knots). 

^ rn3X 

The ratio of ^gA^g = 0.66 was determined by the follcwing procedure 
The single frequency of maximum energy density for each swell train 
recorded at Manterey was read from the f-t curves drawn for each swell 
train as described in Method One. This frequency was converted to for 
fully arisen sea conditions using Figure 8. The geostrophic wind velocity 
at each storm's swell origin point determined by Method Two was obtained 
fron the weather maps. The ratio V^Ag then formed for each storm. 

This set of six values was averaged to obtain the V^Ag ratio used in 
this work. 

2. Results 

Table II presents the results of the procedures followed in Method 
Three to determine the origin and the arrival direction at Monterey of the 
swell trains generated in the five storms studied. 
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3. Time of Arrival of the Peak Eneirgy Swell 

If one assumes that the frequency of the maximum energy swell is 
conserved, the time of arrival, t, of the swell peak at a distant 
station can be estimated using 




vhere is the group velocity of the frequency corponent of maximum 
energy, produced by the surface wind at the peak wind point. 

The results of arrival-time predictions for Monterey obtained in 
this manner for the six swell trains, along with the time of arrival of 
the energy peak as determined frcm spectral analysis of the waves recorded 
at Monterey, are presented in Table III. 

The method tested here is based on an unpublished technique for 
forecasting ocean swell being developed by W. C. Thotpson mder FNWC 
sponsorship. 
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Figure 8: PERIOD AND FREQUENCY OF MAXIMUM ENERGY 

FOR A FULLY ARISEN SEA VERSUS SURFACE 
WIND VELOCITY. (After Pierson and 
Moskowitz , 1964) 
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FREQUENCY OF MAXIMUM ENERGY (HERTZ) 



Table II: RESULTS OF METHOD 
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IV. OCMPARISON OF THE THREE I^ETHODS 



A. OCMPARISCN OF THE THREE METHODS FOR STORM THREE 

As explained in the discussion of Method One, Storm 3 was the only 
storm on viiich all three methods were tested. The peak swell origin 
points as determined by the three methods are shown in Figure 10 by the 
numbered circles. Point (1) marks the approximate intersection of the 
distance arc from Monterey (D = 2330 nautical miles) and the adjusted 
distance arc from San Clemente Island (D = 2574 nautical miles) for an 
origin time of 0300Z 21 November 71 as derived fron wave records at the 
two locations. Based on this intersection the arrival direction of the 
swell at Monterey was 297 degrees true. Points (2) and (3) represent 
origin times of 0330Z 21 Novenrber and 0600Z 21 November 71, respectively, 
and yield swell arrival directions at Monterey of 285. 2°T and 286. 5°T, 
as determined by Methods Two and Three. 

The source origin points for Ifethods Two and Three are quite close 
together (Figure 10) . At first glance the origin point as determined by 
Method One may seem to be poor in relation to the other two; however, 
it should be noted that because of the snail angle of intersection, a 
small change in the source distance determined frcm the f-t curves drawn 
from the spectral analysis of the wave records vail change the north- 
south location of the intersection and the source point by a considerable 
amoxant. Clearly a more suitable baseline between sensors is needed. In 
view of these considerations, the agreement between the three methods 
^plied to this storm is considered to be quite good. 
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Figxire 9; SWELL ORIGIN POINTS FOR STORM 3 AS DETERMINED BY METHODS ONE, 



B. CCMPARISCN OF METHODS TWO AND THREE FOR ALL STORMS 



The ODTtparisons of the swell source distance and origin time, and the 
arrival direction at Monterey, derived using Methods Two and Three for 
the five storms studied are shown in Table IV (obtained fran Tables I 
and II) . The table shows that the peak swell origin time differences 
vary frcm 15 minutes to 6 hours, and the distance differences from the 
origin point to Monterey vary from 37 to 140 nautical miles, with the 
exception of Storm 4A which had a 362 nautical mile difference. The 
differences in arrival direction of the peak swell at Monterey obtained 
fran the two methods varied fran 0.1 to 5.3 degrees. It may be noted 
that Method TWo produced distances greater than those predicted fran 
Rfethod Three in five of the six cases studied. 

Possible reasons for variability in the location of the origin using 
the two methods include: 

1. Subjectivity in the fit of the f-t curve to the spectral energy 
graph. This leads to sane possible variation in the distance and time 
of origin of the observed swell. 

2. Variability in the FNWC weather maps. Lack of data in a storm 
area may lead to an analysis that is unrepresentative of the actual surface 
pressure field present. 

3. Poor estimates of the time of swell origin using Method Three. 

Method Three, because of its sorev^at subjective character, is prone 

to errors in determination of the origin point of the peak swell due to 
errors in estimating the correct time of origin. It is therefore 
appropriate to examine what the effect of a ±12 hour and ±24 hour time 
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error in the choice of time of origin of the peak swell will have on the 
swell arrival direction prediction at a distant coastal station such as 
Monterey. The results of this investigation for the five storms studied 
are shown in Table V. Clearly, an error in the predicted time of origin 
will have little effect on the arrival direction of ocean swell generated 
ky storms at sea that are moving along approximate great-circle paths 
toward the coastal station. However, in the case of a storm such as 
Storm 5, where the path of the storm makes a large angle with the swell 
path to the coastal station, the effects of an error in time of origin 
prediction can be fairly large. It should be noted that from the results 
presented in Table IV, an error in estimation of the swell origin time 
in excess of 12 hours appears unlikely. 
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Table IV: (XMPARISON OF KESULTS FROM METHODS IWO AND 
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Table V: SWELL ARRIVAL DIRECTION AT WDNIEREY RESULTING FROM 

12 AND 24 HOUR ERROR IN ESTIMATING TIME OF ORIGIN 





CO 


VO 




CO 


in 


VD 




• 


• 


• 


• 


• 


• 


CN 




CO 


<j\ 


CN 


CN 


r- 


+ 


o 


o 


CO 


<J\ 


o\ 


iH 




CN 


m 


CN 


CN 


CN 


ro 





CO 




in 


CN 


cr» 




CN 


• 


• 


• 


• 


• 


• 


tH 


in 


CO 


VO 


r- 


o\ 


CO 


+ 


o 


o 


CO 


a\ 


CO 


r— 1 




ro 


ro 


CN 


CN 


CN 


m 





CN 


CN 


o 


in 




• 


• 


• 


• 


• 


• 


VO 


in 


in 


o 


CO 


o 


o 


o 


CO 


o 


CO 


rH 


ro 


ro 


CN 


ro 


CN 


ro 







CN 


1 — ! 


O 


ro 


VO 


CN 


• 


• 


• 


• 


« 


• 


iH 




r- 


VO 


o 


ro 


ro 


1 


o 


o 


CO 


o 


CO 


o 




ro 


ro 


CN 


ro 


CN 


ro 





• 


VO 






r- 






i< 


• 


• 


• 


• 


• 


CN 




r- 


in 


jH 






1 




o 


CO 


o 


CO 








ro 


CN 


ro 


CN 


CN 




42 



1 




LIST OF REFERENCES 



1. Barber, N, F., and F. Ursell, "The Generation and Propagation of 

Ocean Waves and Swell, I. Wave Periods and Velocities," 
Philosophical Transactions, Royal Society of London, v. 240, 
p. 527 to 560, 24 February 1948. 

2. Blackman, R. B., and J. W. Tukey, "The jyfeasureroent of Power Spectra 

from the Point of View of CfcOTnunications Engineering," Dover 
Publications , 1958. 

3. Carstensen, L. P., Fleet Numerical Weather Facility , Technical Note 

No. 29. Seme Effects of Sea Air Terrperature Difference, Latitude 
and Other Factors on Surface Wind Geostrophic Wind Ratio and 
Deflection Angle, p. 1-9, March 1967. 

4. Chang, D., "A Fortran Sub-Routine for the Great-Circle Distance 

Between IWo Points and Bearings at the Points," Naval Research 
laboratory Corrputer Note 32, 1 September 1969. 

5. Lynch, T. J., Long Wave Study of Monterey Bay, Masters Thesis, 

Naval Postgraduate School , Monterey, California, September 1970. 

6. Munk, W. H., G. R. Miller, F. E. Snodgrass and N. F. Barber, 

"Directional Recording of Swell from Distant Storms," Philosophical 
Transactions, Royal Society of London, v. 255, p. 505-584, 

18 April 1963. 

7. Pierson, W. J., and L. Moskowitz, "A Proposed Spectral Form for 

Fully Developed Wind Seas Based on the Similarity Theory of 
S. A. Kitaigorodskii, " Journal of Geophysical Research , v. 69, 
p. 5181-5194, 15 December 1964. 

8. Snodgrass, F. E., G. W. Groves, K. F. Hasselmann, G. R. Miller, 

W. H. Munk, and W. H. Powers, "Propagation of Ocean Swell Across 
the Pacific," Philosophical Transactions, Royal Society of London, 
V. 259, p. 431-497, 4 May 1966. 

9. Thompson, W. C. , "Swell and Storm Characteristics fron Coastal Wave 
Records , " American Society Civil Engineers , Proceedings , Twelfth 
Coastal Engineering Conference, p. 33-52, September 1970. 



43 



INITIAL DISTRIBUTION LIST 



No. Copies 
2 



1. Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 

2. Library; Code 0212 
Naval Postgradiaate School 
Monterey, California 93940 

3. Professor Warren C. Thompson 
Department of Oceanography 
Naval Postgraduate School 
Monterey, California 93940 

4. Lieutenant M. H. Austin, Jr. 
USS Paul Revere (LPA-248) 
c/o FPO San Francisco 96601 

5. Professor J. B. Wickham 
Department of Oceanography 
Naval Postgraduate School 
Monterey, California 93940 

6. Department of Oceanography 
Naval Postgraduate School 
Monterey, California 93940 

7. Lieutenant Commander Charles 
Dep>artment of Oceanography 
Naval Postgraduate School 
Monterey, California 93940 

8. Oceanographer of the Navy 
The Madison Building 

732 N. Washington Street 
Alexandria, Virginia 22314 

9. Dr. Ned A. Ostenso 
Code 480D 

Office of Naval Research 
Arlington, Virginia 22217 

10. Evelyn L. Pruit, Director 

Gejography Programs, Code 414 
Office of Naval Research 
Depjartment of the Navy 
Washington, D. C. 20360 



2 



5 



2 



1 



3 



. Roberts 1 



1 



1 



1 



44 



1 



11. CoriTnanding Officer 
Environinental Prediction Pesearch Facility 
404 Franklin Street 
Monterey, California 93940 

12. Comrnanding Officer 2 

Fleet Numerical Weather Centred 

Monterey, California 93940 

13. Coastal Engineering Research Center 1 

5201 Little Falls Road, N. W. 

Washington, D. C. 20016 

14. Mr. Charles Fisher, Chief 1 

Coastal Engineering Branch 

U. S. Aoty Corps of Engineers 
P. 0. Box 2711 

Los Angeles, California 90053 

15. Ccinmanding Officer 2 

San Francisco District 

U. S. Amy Corps of Engineers 

100 McAllister Street 

San Francisco, California 94102 

Navigation and Shoreline Planning Section 
Library 

16. Coastal Engineering Branch ^ 

Planning Division 

U. S. Amy Engineering Division, South Pacific 

630 Sansome Street 

San Francisco, California 94111 



45 



Unclassified 



S<*cvmty Classification 



DOCUMENT CONTROL DATA - R & D 

[Security classification of title, body of abstract and indexing annotation must be entered when the overall report is c lassified) 



originating activity (Corporate author) 

Naval Postgraduate School 
Monterey, California 93940 



REPOR T SECURI TY CLASSIFICATION 

Unclassified 



2b. GROUP 



3 report title 



Determination of the Deep Water Arrival Direction of Ocean a*;ell at a 
Coastal Station 



4 . DESCRIPTIVE NOTES (Type of report and, inclusive dates) 

Master's Thesis; March 1972 



5- AU THOR (S) fFirs# name, middle initial, last name) 

Marshall Harlan Austin, Jr. 



6 REPOR T D A TE 

March 1972 



7a. TOTAL NO. OF PAGES 



7b. NO. OF REFS 



e«. CONTRACT OR GRANT NO. 



b. PROJEC T NO. 



Oa. ORIGINATOR'S REPORT NUMBER(S) 



9^. OTHER REPORT NOiS) (Any Other numbers that may be aaslgned 
this report) 



10. DISTRIBUTION STATEMENT 

improved for public release; distribution unlimited. 



n. SUPPL EMEN T A R Y NOTES 



12- SPONSORING MILI T ARY ACTIVITY 

Naval Postgraduate School 
Monterey, California 93940 



13. abstract 



Three previously unpublished methods for eitpirically determining deep-water 
swell direction were examined in this work: (1) triangulation using two widely 

separated wave sensors, (2) intersection using weather iraps and a single wave 
sensor, and (3) swell point-source estimation frem weather maps. The primary 
objective of each method was to identify a single point source of the swell 
train produced in an approaching cyclonic storm. Method (1) was not adequately 
tested, but results as applied to swell fresn one storm were favorable. Methods 
(2) and (3) , applied to five selected North Pacific storms, gave close agreement 
on the swell origin time (15 minutes to 6 hours) and source point location 
(27 to 362 nautical miles) , and on swell arrival direction at Monterey, California 
(0.1 to 5.3 degrees). Methods (1) and (2) give deep-water directions for swell 
trains that have already arrived, while Method (3) is a prediction method. 



DD, 



fo«. 1470 

NOV 05 I / hj 

S/N 01 01 -807-681 1 



(PAGE 1) 



Unclassified 



46 



Security ClassificBtion 



A-31408 




J 



◄ 












m 




Unclassified 

Security Classification 



KEY wo R DS 



ROLE W T 



ROLE W T 



ROLE w T 



Swell Propagation 
Frequency-Time Graph 
Swell Arrival Direction Prediction 
Great-Circle Path 
Great-Circle Distance 
Great-Circle Azimuths 



FORM 

1 NOV 

101 -807-6821 



..1473 ( BACK ) 



Unclassified 



47 



Security Classification 



A- 3 I 409 



f 



s 

I 

[ 

i 

j 

1 

( 

I 



'Thes i s 
A9715 



TV 

A 



Dupl i cate 

133999 

Austin 

Determination of the 
deep water arrival di- 
rection of ocean swell 
at a coastal station. 



Thesis 
A9715 
c .1 



~ 133999 

Austin 

Determination of the 
deep water arrival di- 
rection of ocean swell 
at a coastal station* 



